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ABSTRACT 

It is known that vacuum polarization can modify the photon propagation modes in the atmospheric 
plasma of a strongly magnetized neutron star. A resonance occurs when the effect of vacuum polariza- 
tion on the photon modes balances that of the plasma. We show that a photon (with energy E ^ a. few 
keV) propagating outward in the atmosphere can convert from one polarization mode into another as it 
traverses the resonant density, p^^s — Y~^ri~^{B /lO^'^ G)'^{E/1 keV)^ g cm~'^, where Yg is the electron 
fraction and 77 1 is a slowly-varying function of the magnetic field B. The physics of this mode conver- 
sion is analogous to the Mikheyev-Smirnov-Wolfenstein mechanism for neutrino oscillation. Because the 
two photon modes have vastly different opacities in the atmosphere, this vacuum-induced mode conver- 
sion can significantly affect radiative transport and surface emission from strongly magnetized neutron 
stars. 

Subject headings: magnetic fields - radiative transfer - stars: neutron 



1. INTRODUCTION 

Soft gamma-ray repeaters (SGRs) and anomalous X-ray 
pulsars (AXPs) are enigmatic objects. It has been sug- 
gested that these are neutron stars (NSs) endowed with 
superstrong magnetic fields, B ^ 10^'* G (see Thompson & 
Duncan 1996). Support for the "magnetar" hypothesis has 
come from several recent observations (e.g., Kouveliotou 
et al. 1998; HuUeman et al. 2000; Kaplan et al 2001; Kaspi 
et al. 2001; see Hurley 2000 for SCR review; Mereghetti 
2000 for AXP review). Thermal radiation, attributable to 
surface emission at T^s = (3 — 7) x 10^ K, has already been 
detected in four of the five AXPs and in SCR 1900-^14 (see 
Mereghetti 2000; Perna et al. 2001). These surface X-ray 
emissions can potentially reveal much about the physical 
conditions and true nature of SGRs and AXPs. Theoreti- 
cal studies on the atmospheres and radiation spectra from 
magnetars are still in the early stages (see Ho & Lai 2001a 
and references therein). 

In this paper, we are concerned with the effect of vacuum 
polarization on the radiative transport in the atmospheric 
plasma of strongly magnetized NSs. It is well known that 
in strong magnetic fields, vacuum polarization due to vir- 
tual e^e~ pairs can modify the dielectric property of the 
plasma and the polarization of photon modes, thereby al- 
tering the radiative opacities (e.g., Adler 1971, Tsai & 
Erber 1975; Meszaros & Ventura 1979; Pavlov & Gnedin 
1984; Heyl & Hernquist 1997; see Meszaros 1992 for re- 
view). Of particular interest is the "vacuum resonance" 
phenomenon, which occurs when the vacuum and plasma 
effects on the linear polarization of the modes cancel each 
other, giving rise to a "resonant" feature in the opacities 
(e.g., Pavlov & Shibanov 1979; Ventura et al. 1979). Re- 
cent works have begun to address this opacity effect of vac- 
uum polarization on the radiation spectra from magnetars 
(Bezchastnov et al. 1996; Bulik & Miller 1997; Ozel 2001). 
At a given density p (in g cm^"^), the vacuum resonance is 

^ This adiabatic mode conversion is not to be confused with the physic 
(see Meszaros 1992), which has been included in previous works (e.g., 



located at photon energy Ey — l.O {Ye pY^^B^^r] keV, 
where Yg = Z/A is the electron fraction and 77 ~ 1 
is a slowly-varying function of the magnetic field B [see 
eq. (11)]. Because Ey depends on p (which spans a wide 
range in the atmosphere), one expects that any spectral 
feature associated with Ey will be significantly spread out 
in the spectra of magnetars (Ho & Lai 2001b). 

Here we study a new physical effect overlooked by all 
previous works: a photon propagating down the density 
gradient in the NS atmosphere can change from one po- 
larization mode to another at the vacuum resonance. This 
conversion is particularly effective for photons with E' <; a 
few keV. The physics of such mode conversion is analogous 
to the Mikheyev-Smirnov-Wolfenstein (MSW) neutrino os- 
cillation which is thought to take place in the Sun (e.g., 
Bahcall 1989; Haxton 1995)^ We show that the mode 
conversion can dramatically affect radiative transport in 
magnetar atmospheres. 

2. PHOTON MODES IN A MAGNETIZED ELECTRON-ION 
PLASMA INCLUDING VACUUM POLARIZATION 

In a magnetized plasma -I- vacuum medium, the dielec- 
tric tensor e and inverse permeability tensor p, can be writ- 
ten as the sums of the usual plasma terms and the vacuum 
correction, i.e., e = e^^^ -I- Ae'^^ p, = I + A/i'^"^ (where 
I is the unit tensor). The plasma (electrons and ions) di- 
electric tensor e'^^^ is given in, e.g. Shafranov (1967). The 
vacuum corrections can be written as 

Ae(^) = (a-l)I+qBB, A/i^"^ = (a-l)I+mBB, (1) 

where B is the unit vector along the magnetic field B, and 
a, q, and m are functions of B. For B <C Bq = mlc^ /eh = 
4AU X 10^3 Q 



a ~ 1 — 25, q — 75, m = —45, with 5 — 




different (and less drastic) effect of mode switching due to scattering 
lov et al. 1995; Ho & Lai 2001a; Ozel 2001). 
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where a = 1/137 and b = B/Bq (Adler 1971). For 
B ^ Bq, using the expansion in Heyl & Hernquist (1997) 
(see also Tsai & Erber 1975), we have 
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An electromagnetic (EM) wave with a given frequency w 
satisfies the wave equation 

, ,2 



V X (jii • V X E) = ^ € • E. 



(6) 



For normal modes propagating along the 2:-axis, with 
E oc e*'^*^, this reduces to 



-k\_z X [/i • (i X E±) 



1^ ^ 
— e-E±, 



(7) 



where the subscripts "±" specify the two modes ("plus- 
mode" and "minus- mode" ) . Using the procedure of 
Meszaros & Ventura (1979; who neglected ions and con- 
sidered the 6^1 limit), the mode eigenvector E± can be 
expressed explicitly: In the xyz coordinates with k along 
the z-axis and B in the x-z plane (such that A: x B = 
smOsy, where is the angle between k and B), we write 
E± = E±T + E±zZ, and the transverse part is 

E±T^ , i,,,,J zj^±,l), (8) 



(1 + ii:j)V2 



where 



K±=p± + r. 



(9) 



with r = 1 + {m/a) sin^ 9b — 1, and (assuming v <^1) 



1/2 



(1 - Uj)sin^ 6b 



2 cos 6Ib 
Here Ue = (wse/w)^ 



1 - 



q + m 



(1 



(10) 



and the electron cyclotron frequency lobc, the ion cy- 
clotron frequency ubi, and the electron plasma frequency 



are given by hujBe 



hujBi 



-- hZeB/{AmpC) 



heB/{mec) = 1158^14 keV, 
= 0.6305 Si4 keV, hwp = 



28.71 (Yepi)^/^ eV. We shall be in- 
terested in photon energies such that Ug ^ 1 and v ^ 1. 

The standard way of classifying photon modes in a mag- 
netized plasma (e.g., Meszaros 1992) relies on the value of 
A' I = \Ex/Ey\: the extraordinary mode (X-mode) has 
I AT I <C 1, and its E is mostly perpendicular to the fc-B 
plane; the ordinary mode (0-mode) has \K\ ^ 1, and 
is polarized along the fc-B plane. The advantage of such 
classification is that the X-mode and 0-mode interact very 
differently with matter: the 0-mode opacity is largely un- 
affected by the magnetic field, while the X-mode opacity 
is significantly reduced (by a factor of order oj^/ojg^) from 
the zero-field value. However, this classification becomes 
ambiguous when <; 1 (\K\ ^ 1). Obviously, f3 — 
for u) = LOBi- But even for general energies {E ^ hcuBi), 
a photon traveling in an inhomogeneous medium encoun- 
ters /? = when the vacuum resonance {v = q+m) occurs, 
namely at the resonance density 

p,es = 0.964 y-i(Si4i;i)'r?-2 g cm-^ (11) 



where Ei = E/ (1 keV) and rj = [36 /{q + m)]^/^ is a slow- 
varying function of 6 [ry = 1 for 6 4; 1 and r] — > (6/5) = 
0.673 bH"^ for 6 > 1; varies from 0.99 at Bu = 1 to 6.7 
at Bu = 100]. 

Figure 1 shows the values of K± and the refractive in- 
dices n± = ck±/uj for the plus and minus-modes near the 
resonance density. The plus-mode (minus-mode) mani- 
fests as O-modc; (X-mode) at high densities but becomes 
X-mode (0-mode) at low densities. A natural question 
arises: how does the polarization state evolve as a photon 
traverses from the high density region through p^cs to the 
low density region? Before we proceed, it is convenient to 
introduce the "mixing" angle 9m via tan^^ = so 
that 

tan2em = P~\ (12) 

where we have used r — 1 <C 1. The eigenvectors of 

the modes are E+t = (* cos6'm,sin6'm) and E_r = 
{—ism 9m, cos 9m)- Clearly, at the resonance 9 m = 45°, 
the X-mode and 0-mode are maximally "mixed" . 

3. MODE CONVERSION IN AN INHOMOGENEOUS PLASMA 

If the density variation is sufficiently gentle, an 0-mode 
photon created at high densities will remain on the 
trajectory (Fig. 1) as it travels outward and will adiabat- 
ically convert into the X-mode after traversing the reso- 
nance density. This is analogous to the MSW mechanism 
for neutrino oscillation. To quantify this, we write the 
general, polarized EM wave with frequency u) traveling in 
the ^-direction as 

F,{z) = A+-E+e''^+ + A_E_e*"^- , (13) 

where (f>±{z) = J^' k± dz. We shall adopt the WKB (or ge- 
ometric optics) approximation, in which the A± and E± 
vary on a length scale much longer than the photon wave- 
length. Substituting eq. (13) into eq. (6), using eq. (7) 
and the approximations a — 1 ^ 1 and m ^ 1, and then 
multiplying the resulting equation by E^^ and E^^, re- 
spectively, we obtain the amplitude evolution equations 
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A+ 



k-m*^ 2k- 

(14) 

where ' stands for d/dz and = ^_|_ — Neglecting the 

off-diagonal terms, we find A±{z) oc \k±{z)\~^/'^ , which is 



A±e"f'±, 



we can 



a standard WKB result. Defining Aj 
cast eq. (14) into a form analogous to that of quantum 

mechanics: 

a; \ / -Ak/2 i9'm \( A+ 
-\ -i9L Ak/2 j I A- 



(15) 



where Afc 



A:_ , and we have subtracted a nonessen- 



tial unity matrix and used k-^ ~ fc_ and |fcj_/fc±| <C \k±\. 



Clearly, when \9!^ 



< |Afc/2|, or 

(n+ 



29L 



>1, 



(16) 



the polarization vector will evolve adiabatically (e.g., a 
photon in the plus-mode will remain in the plus- mode). 
Using eqs. (12) and (10), we find 



isin^ 29„ 



1/2 



(1 



u, sm 



2 cos 9b 



^ U . (17) 
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The difference in refractive indices of the two modes is 

vcos 6b 



n+ — n_ ~ — - 



1/2 



(1 — Ui) sin2^„ 



(18) 



Clearly, the adiabaticity condition (16) is most easily vio- 
lated at the resonance {9m = 45°). Evaluating eq. (16) at 
P = Pres, we obtain 

7,es = {E/E^f, (19) 

1/3 



= 2.52 (77 tan 



2/3 



1 



|l/3 



1 cm 



keV, 



(20) 

where Hp = \dz/dhip\ is the density scale height along the 
ray (evaluated at p = pres)- For an ionized hydrogen at- 
mosphere, Hp ~ 2kT/{mpg cos 6) = 1.65 Tg/ (514 cos ^^) cm, 
where T = lO^Tf, K is the temperature, g = 10^"'.gi4 cm s~^ 
is the gravitational acceleration, and 9 is the angle be- 
tween the ray and the surface normal. Thus, for E ^ 
2Ead, resonant conversion between X-mode and 0-mode 
is essentially complete; for E ^ 0.5£^ad, a photon will 
"jump" across the adiabatic curves (Fig. 1), and an X- 
mode (0-mode) will remain an X-mode (0-mode) in pass- 
ing through the resonance. In general, the jump probabil- 
ity can be calculated with the Landau- Zener formula 

p. _ p-7r7rcs/2 (2-1) 

^ jump — ^ v^^/ 

(see Haxton 1995). The width of the resonance region can 
be estimated by considering = 1 as defining the edge 
of the resonance; this gives 

(A^).es ^ 1.7 X 10-3|1 - u,\-'^^ Hp, (22) 

sm^ 9b \-di4 / 

which is much smaller than the photon mean-free path. 

4. IMPLICATIONS FOR SURFACE EMISSION FROM 
MAGNETARS 

To understand qualitatively the effects of vacuum po- 
larization (and mode conversion, in particular) on the ra- 
diation spectra from magnctar atmospheres, wc begin by 
estimating the location of the decoupling layer (at which 
the optical depth is 2/3) for photons of different energies 
and polarizations. For simplicity, we neglect scattering 
and consider fully-ionized H atmospheres. 

First, let us turn off the vacimm polarization and pro- 
ton effects. The free-free absorption opacity can be written 

as K = KoC, where kq = 9.3piTg ^^^E^'^G cm^g^^ (with 
G — 1 — e"^/*^"^) is the zero-field opacity and ^ ~ 1 for the 
0-mode and ^ ^ for the X-mode (we set the Gaunt 
factor to unity). Hydrostatic equilibrium yields the col- 
umn density y ~ 0.83 piTq g cm~^ (for 514 = 2). The 



decoupling densities of the two modes are then Pq 



(nv) 



(nv) 



1/2 r 

po Ue I see 



atmosphere, and thus they are the main carriers of the 
X-ray flux. 

Next we include vacuum polarization and protons, but 
turn off mode conversion (e.g., X-mode remains X-mode 
across the vacuum resonance). The decoupling density 

po for the 0-mode is still pQ^^ ■ For the X-mode, when 

(nv) \ 
PX <' Pres, or 

(23) 



E^E,2^ 63 76 '/'(Bi4/ry2)-2/3 keV, 



the photons created at p^'^' will encounter the resonance, 
near which the X-mode opacity is greatly enhanced; thus 

the decoupling density p^'^'' ["(nc)" stands for "no conver- 
sion"] will be smaller than p^^'' (see Fig. 2; for E <; Ec2, 

p'^^^ is close but not equal to p^^-* because the opacity 
is modified slightly by the vacuum effect even away from 
the resonance). In fact, the optical depth across resonance 



region (say 0.9 < p/pres < 1-1)^ is a factor 



,1/2 



larger 



0.42Te-^/"i;f/'G-i/2 g cm'^ and pl^ 
Fig. 2; the superscript "(nv)" stands for "no vacuum"; 
note that we have treated the temperature as a constant; 
this a good approximation when estimating the decoupling 
density since T varies by only a factor of a few while p 
varies by many orders of magnitude above the decoupling 
layer]. The X-mode photons emerge from deeper in the 

The X-modc opacity is k = kqC, witli g = sin^ 9^/(2/3)'^ + uj^ for | 
|a;| = |Ap/pres| ^ 0.1, we have j3 ~ u\^^x. Thus £, ~ l/{u^x'^) for Ue ^^'^ 

— 1/2 

the ^ 2ue region contributes most to Ar. 



than that of the non-resonant region (p <, 0.9pres)- Thus, 

as E decreases below Ec2, we expect p^^'^^ to follow closely 
Pres) until E drops below another critical energy Eci (see 
Fig. 2), which is set by Ar ~ 2/3, or 

^ei ^ 6 T-^'^B-^'^-" keV (24) 
(the numerical prefactor 6 is for Ob — 45°). Below Ed, we 

find that pl^'^' becomes close to p^^'' again, except for the 
proton cyclotron feature at Ebi = 0.63 Bn keV (studied 
previously by Ho & Lai 2001a and Zane et al. 2001). 

Now consider the effect of mode conversion at the vac- 
uum resonance. For E > 1.3£^ad) adiabatic mode con- 
version is nearly complete (Pjump < 3%). The 0-mode 
photons traveling from high densities through prcs are con- 
verted to X-mode photons, which then freely stream out 
of the atmosphere. Thus for 1.3 E^ ^ E ^ Ec2, the ef- 
fective decoupling density for X-mode is px = Pres- For 
E Ec2- vacuum resonance occurs inside both the X 

and 0-mode decoupling layers, and therefore px — p'x^^ • 
For E ^ Eg^, mode conversion is ineffective, and thus 
Px = Px'^^- Around E'ad, the X-mode photons are emit- 
ted from both p'^'^^ (with probability Jjump) and Pres [with 
probability (1 - Pjump)]- 

To translate Fig. 2 into emergent spectra from NS 
atmospheres, we need to know the temperature profile 
T(p). Obviously, T(p) can only be determined from 
self-consistent atmosphere modeling, and different pho- 
ton decoupling behaviors will give rise to different tem- 
perature profiles. This is beyond the scope of this pa- 
per. Here, for illustrative purposes, we consider a fixed 
temperature profile, obtained from the B = 10^^ G, 
Tcs ~ 5 X 10^ K, hydrogen atmosphere model (which in- 
cludes the proton effect but not vacuum polarization) of 
Ho & Lai (2001a; see Fig. 4 in that paper); this profile is 
strongly non-Eddington, showing a mild plateau around 
p = 2 — 20 g cm~^ {Tq varies from 3.5 to 3.7), as a result 
of radiative transport by two modes with vastly different 
opacities. Given this T(p), we calculate the emergent spec- 
tral flux by ~ (7r/2) {B,[T(px)] + B.[T(po)]}, where 

/3| > 1 and ^ = siii-(9B/2 for |/3| < 1 (sec Ho & Lai 2001a). For 



■C |a;| ^ 0.1 and $ ~ 1 for |x| < u, 



-1/2 



(for typical Ob's). Note that 
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B^(T) is the Planck function and px, po are given as in 
Fig. 2. In the case of partial mode conversion, we replace 

B,[T{px)] by (l-Pjump)S,[r(pres)]+PjumpB,[T(p^"))]. 

Figure 3 shows the results for the different cases discussed 
in the previous paragraphs. We see that the vacuum polar- 
ization reduces the flux in the energy band E^d ^ E ^ Ec2- 
Part of the depletion in Fig. 3 is caused by our assumption 
of fixed temperature profile (e.g., the heavy-solid curve and 
the dashed curve have different total fluxes), but Fig. 3 is 
indicative of the potentially important effects of vacuum 
polarization (and mode conversion) on the radiation spec- 
tra from magnetars. 



For Pq^^ ^ pres, i-e., i?i4 ^ 0.66Ye^^ E^ the vacuum 
resonance occurs outside the O and X-mode decoupling 
layers, thus the total spectral flux {Fx + Fq) is largely 
unaffected by the vacuum resonance, although the polar- 
ization will still be affected. 

Previous studies (e.g., Ozel 2001; see footnote 1) ne- 
glected adiabatic mode conversion as discussed in this pa- 
per. We expect that including this effect in the radiative 
transport would produce a qualitatively different spectrum 
for magnetars (see Ho & Lai 2001b for details). 

This work is supported in part by NASA grant NAG 
5-8484 and an Alfred P. Sloan fellowship to D.L. 
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Fig. 1. — The polarization parameters (upper panel) and refractive indices (lower panel) of photon modes as functions of density near the 
vacuum resonance for B14 = 1, E = 1 keV, Ye = 1, and 6b = 45°. 
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Fig. 2. — The decoupling densities (at which r = 2/3) of different polarization modes as a function of photon energy for B14 = 1, 6b = 45°, 
and T ~ 5 X 10^ K. The solid lines show the case when vacuum polarization is included but mode conversion is turned off; the short-dashed 
line shows the case when the vacuum polarization and proton effects are turned off (note po is unaffected by vacuum effect); the long-dashed 
line shows the vacuum resonance density. The thick vertical line corresponds to the critical energy for adiabatic mode conversion. 
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Fig. 3. — Spectral fluxes from NS atmospheres at B14 = 1. All models assume the same temperature profiles based on the Tefi = 5 x lO'' K 
model of Ho & Lai (2001a). Note that the results are deduced from Fig. 2 (see text) and are qualitative; they serve to illustrate the effects of 
vacuum polarization and mode conversion. 



